Heat pumping with optically driven excitons 
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We present a theoretical study showing that an optically driven excitonic two-level system in a solid state 
environment can act as a heat pump by means of repeated phonon emission or absorption events. We derive a 
master equation for the combined phonon bath and two-level system dynamics and analyze the direction and rate 
of energy transfer as a function of the externally accessible driving parameters in the coherent control regime. 
We discover that if the driving laser is detuned from the exciton transition, cooling the phonon environment 
becomes possible. 



In semiconductor quantum dots (QDs) the ground state 
and the state containing a trapped electron-hole pair (exciton) 
form a two level system (2LS) which is a popular implemen- 
tation of a qubit. Unlike their atomic counterparts, such ex- 
citonic qubits are inextricably coupled to the lattice dynamics 
of the surrounding material fl}|5). Longitudinal optical (LO) 
phonons give rise to excitonic polarons 1 3 ] and are known to 
play an important role for ultrafast excitation [6, 7]. However, 
in the much slower coherent control regime that is so inter- 
esting for quantum information processing, optical phonons 
only contribute negligibly to the dephasing, and the coupling 
to longitudinal acoustic (LA) phonons via the deformation po- 
tential is predicted to be dominant at low temperatures |7]|8). 
In recent experimental studies, the influence of LA phonons 
on the coherent QD dynamics has been measured by the tun- 
nel charge through the QD (proportional to the excitonic pop- 
ulation) following optical excitation [9], and also through the 
direct coupling of the QD's dipole to optical modes in mi- 
crocavities ffTOlfTTI . These studies confirm that a 2LS model 
with a perturbative treatment of the LA phonon interaction 
provides an excellent description in this regime. 

Previous theoretical studies of the exciton-phonon interac- 
tion in the Rabi regime have fully discarded all information 
about the state of the phonon bath l6ii8l [12UI41 . Here we de- 
velop a technique for tracking the number of excitations in the 
environment, allowing us to analyze the net rate of absorbed or 
released bath energy, and we show that a continuously driven 
excitonic qubit constitutes a controllable two-way heat pump. 
Exploiting this cooling effect would help with gaining further 
experimental insight into the exciton-phonon coupling, crucial 
for managing decoherence of semiconductor charge qubits, as 
well as providing an easy preparatory qubit initialization step 
for quantum information processing. 

Model - We consider a self-assembled QD (such as In- 
GaAs encased in GaAs substrate) illuminated by a laser beam 
with frequency oji, which is nearly resonant with the crystal 
ground state to exciton transition. In a frame rotating with the 
laser frequency and within the rotating wave approximation 
(RWA), the system is governed by the Hamiltonian (h = 1) 
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FIG. 1: (color online) The driven exciton acts as a heat pump be- 
tween the phonon and the electromagnetic environment. If the driv- 
ing is detuned from the exciton transition, the laser-dressed eigen- 
states (see Ref. 1 14j) are composed of different amounts of the 
ground and excited state. Cooling is possible when spontaneous 
emission repeatedly takes the system into its lower eigenstate, from 
where it can only absorb phonons. 



coupling the two basis states, and where a x = \g){e\ + \e)(g\ 
and a z — |e)(e| — \g)(g\ are the Pauli pseudo-spin operators 
defined with respect to the QD ground state \g) and single- 
exciton state |e). 

Let us further assume that the excitonic qubit is coupled to 
a bath of phonons, resulting in the total Hamiltonian 



H — Hq + H r + Hi 



(2) 



where Hb = 53 q w q a q a q is the free Hamiltonian of the 
phonon modes with a q and a q being the creation and anni- 
hilation operators of a phonon in mode q with frequency w q . 
The exciton-phonon interaction term is generically given by 



Hi = cr z ^5 q (a q + o q ), 



(3) 



where the <7 q are coupling constants. Moving into the inter- 
action picture with respect to H$ and Hb we obtain for the 
transformed interaction Hamiltonian 



H s = A/2a z + n/2a x , 



(1) 



where the detuning A describes the energy difference between 
the basis states in the rotating frame, is the Rabi frequency, 



where a z (t) denotes the transformed a z operator. 
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For monitoring phonon excitations in the QD's solid state 
environment, we adapt a technique from the literature on ob- 
serving single electron charge transport |[T5l4T7ll by perform- 
ing a gauge transformation which adds phase markers to the 
interaction Hamiltonian. Crucially, these do not alter the sys- 
tem's dynamics. We define the phonon-number specific den- 
sity matrix, i.e. an object that contains information about both 
the qubit dynamics and the number of phonons in the environ- 
ment, 



Pn 



tr E (P m g), 



(5) 



where trg denotes the trace over the phonon modes, g is the 
density matrix of the combined phonon-qubit system and P ro 
is a projection operator, 
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with the total phonon number operator N = ^ q u q u q 
the gauge transformation operator £\. The projection op- 
erator P m projects out all parts of the wavefunction with a 
phonon number different from m (for a detailed discussion 
see Ref. |fl6l ). The probability of having emitted rn phonons 
into the environment is then obtained by the matrix trace of 
the phonon-number specific density matrix 



Pn 



tr(p„ 



(6) 



It is convenient to introduce the Fourier transformed density 
matrix as 



E 



imX 



with the definition q\ = £x/2Q^\/ 2 - This particular choice 
corresponds to an initial state with a definite phonon number 
as in Ref. 1151 . It can be shown that the full density matrix g\ 
in the interaction picture obeys the von Neumann equation, 

Q\{t) = -% (ilx{t)Q\{t) - Qx(t)H- X (t)) = Cx(t)g x (t), 

defining the super-operator C\(t), where the interaction 
Hamiltonian H\ = £\/ 2 Hi£^_ x / 2 has acquired phase markers 
on the phonon creation and annihilation operators: 



H 



x 



q 



( 7 ) 



written in a shorter notation as H\ — d z (t)B\(t). The 
phase factor exp(— iA/2) then keeps track of the creation of 
phonons, while exp(zA/2) tracks annihilation processes. 

We proceed along the standard path of deriving a master 
equation (ME) |18|, and obtain an integro-differential equa- 
tion for the reduced density matrix of the system p\, 

t 

px(t) =tr E / ds£\(t)£\(s)g\. 



After the Born-Markov approximation [18] the resulting 
Markovian ME reads 



Px(t) 



poo 

/ ds^GiM^Mt-^Pxit), 



where we have defined the super operators Si(t)X = 
a z (t)X, S2(t)X = X<r z (t) and the environment cor- 
relation functions G^(s) = (P A (s)B A (0)}, G x 2 (s) = 
-(B-x(0)B x (s)), G^(s) = -(B_ x (s)Bx(0)) and 
G x 2 (s) = (P_ A (0)P_ A (s)). To simplify the further alge- 
braic evaluation, we express the operator a z (t) in terms of 
system eigenoperators, yielding 



E 

ue{o,A} 



'e~ iut P^ 



Jut pf 



Pi) 



(8) 



where A = \/il 2 + A 2 now denotes the spacing between 
the system eigenstates |— ) = cos 6%) — sin0|e) and |+) = 
cos9\e) +sin 0\g) of Hamiltonian ([T]i with 29 = arctanf2/A. 
In this basis, P = cos 26(\-)(-\ - |+)(+|)/2 and P A = 
sin 26>|— ) (+|. 

By introducing the Hermitian operator V = 2Po + 
e -»Atp A _j_ e iAtpt an( j abbreviating Q = e~ lAt P\, we obtain 
the following interaction picture ME after some straightfor- 
ward algebra 

Px = T x (V iA (QrxP^ +f P x&) - V^Qpx - px&V 
+ T t (e lX (v^pxQ+ &PxV) - V&PX - PxQVi 



where the rates are given by = J (A) (n(A) + 1) /2 and 
= J(A)n(A)/2 when the phonon bath pe is in a thermal 
state ED. Further, 



J(w) = 2tt V |g q | <5(w-w q ) 



q 



(9) 



is the spectral density of phonon modes and n(u>) — 
(exp(/3w) — denotes the thermal occupancy of a phonon 
mode with frequency uj. 

Upon going back to the Schrodinger picture, the operators 

V and Q lose their time-dependent phase factors, turning into 

V = a z and Q = P\, respectively. We can return to the 
number representation with the transformation 
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This yields set of coupled differential equations for the evo- 
lution of phonon-specific density matrices. After a RWA (im- 
plemented by setting V = C, expected to be valid whenever 
A > J (A) QUE]]), we finally obtain a ME in diagonal Lind- 
blad form 

pm = -i[Hs(t),p m ] 

+ T± (2P APm+1 P A - PlP A p m - PmPlPh) 



r t (2P A p m _lP A - P A P{ Pm - PmP^Pl 



■(10) 
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This equation describes the joint phonon-qubit dynamics in 
the case where the phonon environment is only weakly dis- 
turbed from thermal equilibrium by the excitonic qubit. 
Phonon-assisted transitions - We now apply Eq. 



10 1 to an 



excitonic qubit that is optically driven by a pulse of constant 
intensity fl. The initial state at t = is the system ground state 
with zero excitations in the environment: p n (0) = S n fi\g)(g\. 
For simplicity, we define the spectral density phenomenolog- 
ically: J(u>) = aui 3 exp (— cj 2 /u; 2 ), where a describes the 
effective electron-phonon coupling strength and uj c is the high 
frequency phonon cut-off. For relatively weak driving with a 
peak Rabi frequency A well below both the electron and the 
hole cut-off, we can neglect the exponential cut-off term alto- 
gether. Setting a = 1/4 ps 2 yields a coupling strength that 
is consistent with the magnitude of the GaAs deformation po- 
tential reported in the literature 



The structure of Eq. ( 10 1 permits the emission or absorption 
of no more than a single phonon: The system is initialised in 
the state \g)(g\ = (|-) + |+))«-| + (+|)/2, i.e. in an equal 
superposition of system eigenstates. The Lindblad operator 
P\ induces a transition from |+) to |— ) while Pi raises pop- 
ulation from |— ) to |+). Once a decay process has happened, 
we find the system in the |— ) state, meaning it cannot decay 
again. 

Provided the excitation is sufficiently long, the population 
ratios thus tend to a Boltzmann distribution, as is obvious from 
the phonon emission rate proportional to n(A) + 1 and the 
absorption rate proportional to n(A), 
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So far, the only perturbation to the system has been caused 
by the coupling to the phonon bath, resulting in single phonon 
emission and absorption. Realistically, other dissipative pro- 
cesses will be present in any physical systems. We shall in- 
clude these using additional phenomenological noise opera- 
tors, such as pure dephasing and radiative decay of the exci- 
ton. We model these processes with an additional Lindblad 



dissipator [18] on the right-hand side of Eq. ( 1 1, 



Vp = T(LplS-±(tfLp- 



(11) 



with respective Lindblad operators L — <j z and L = er_, and 
where T is the dephasing or decay rate. These operators do not 
preserve the system's eigenstates; consequently, under their 
action, phonon-assisted transitions become possible in any of 
the p n subspaces. This leads to a dynamic equilibrium, where 
phonon-assisted transitions keep occurring after the transient 
(coherent) evolution of p = J2 n P n rias subsided. Our theoret- 
ical model is well-suited for illustrating this behaviour: Fig. [2] 
presents the p n distribution at different points of time for the 
case of optical decay. 

The pure dephasing a z operator randomises the phase be- 
tween \g) and |e), thus balancing the population of |— ) and 
|+). Consequently, once the steady state has been reached, 
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FIG. 2: (color online) The distribution of p„ at different points of 
time in the presence of radiative decay. The constant excitation pulse 
uses Q = 1 ps -1 (~ 2/3 meV) and a detuning A as indicated in 
the plots (in ps -1 ). The decay rate is fixed at T — 0.1 ps -1 . The 
different colors/symbols correspond to time as follows: red squares 
= 1.2, green circles = 10, blue diamonds = 40, gray triangles = 70, all 
in units of full Rabi cycles, 2n/\ / 'Q 2 + A 2 . 



phonon emission always occurs with a faster rate then absorp- 
tion, and the distribution is therefore shifted in the direction of 
increasing n, meaning the average number of emitted phonons 
increases. On the other hand, the distribution can move in ei- 
ther direction for optical decay from |e) to |<?). For A = 0, 
\g) consists of an equal superposition of |— ) and |+), while 
it contains a larger |— ) component for |A| > 0. Under these 
latter circumstances, it is possible for phonon absorption to 
permanently dominate over emission, shifting the distribution 
in the direction of decreasing n, as shown in the right panel of 
Fig. [2] In this case, thermal energy is removed from the QD's 
bulk surroundings and released into the wider environment by 
spontaneous photon emission (depicted in Fig. [TJ. 

Heat transfer rate - To quantify this observation, we con- 
sider only the radiative decay operator of Eq. ( fTT| in the fol- 
lowing. We proceed by calculating the rate of phonon emis- 
sion or absorption, which is given by ifTTl : 



(n(t)) 



d 



m 

2tr(r^ep(*)et-r t eV( t ) e ) ) ( i 2 ) 



A=0 



where pit) is obtained by integrating Eq. 
Lindblad dissipator (111 with L = er_ 



( |10| inclusive of the 
and disregarding the 
indices m of p m . 

Fig. [3]presents the steady-state value (h) ss of Eq. (12i as 
a function of A. As expected, net phonon absorption is only 
possible at finite temperature for off-resonant excitation. It is 
a natural question to ask at which rate energy is transferred to 
or from the surroundings of the system. Considering the quan- 
tity E — y/il 2 + A 2 (h) ss shows that the number of absorbed 
phonons decreases for a larger detuning, yet their energy is 
greater, shifting the A that achieves optimal heat transfer. The 
inset of Fig. [3] shows that the process is limited by the ra- 
diative decay time T of the system. A saturation only occurs 
when the spontaneous emission rate becomes faster than that 
of phonon-mediated transitions, but this regime would require 
an unrealistic optical lifetime of the order of a picosecond or 
less. 
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FIG. 3: (color online) The rate of phonon-induced transitions {h) ss 
(red) and the rate of energy transfer E (blue), both as a function of the 
detuning for Q = 1 ps -1 . The rates shown here are the steady state 
values of Eq. 1 12 1 for the decay rate T = 0.1 ps -1 . A positive sign 



indicates net phonon emission whereas a negative corresponds to net 
phonon absorption. The inset shows {n) ss for fixed A = 1 ps -1 as 
a function of T at T = 10 K. 



We proceed by estimating the achievable cooling rate for 
realistic parameters in SI units. Using f2 = A = 1 ps -1 
corresponds to a phonon energy of A = 1.49 x 10~ 22 J. 
With approximately 0.02 absorbed phonons per picosecond 
(at T = 20 K), we obtain a theoretical energy transfer rate 
of roughly 3 x 10~ 12 J/s. Neglecting heating effects, this 
achieves a temperature reduction of the order of one Kelvin 
per second for a micrometer cube of GaAs [22]. Of course, 
any real sample will also be subject to heating processed, 
e.g. by the laser illumination and by thermal contact to its sur- 
roundings. While it is difficult to estimate the precise rate 
of heating, it could plausibly exceed the cooling rate in bulk 
samples. However, by incorporating the QD on a nanopillar 
or a cantilever, the thermal coupling to the bulk could be re- 
duced to make net cooling possible ||23ll . An experiment to 
show the cooling effect could proceed in the following way: 
insert a cooling cycle before the standard sequence for observ- 
ing Rabi oscillations and study the influence of the length of 
the cooling cycle on coherence times. If cooling is successful 
one would expect increased coherence times. 

Summary - We have shown that a single excitonic QD, an 
experimentally well studied system, can act as a heat pump. 
Adapting a counting statistic approach from the context of 
charge transfer, we have discussed how energy can be re- 
moved from the QD's environment by means of repeated 
phonon absorption in conjunction with spontaneous photon 
emission. This effect does not rely on the structure of the spec- 
tral density, making our analysis applicable to similar systems 
with a different coupling mechanism to a bath of harmonic 
oscillators. This opens up the possibility of experimental in- 
vestigation of quantum heat pumps as well as the prospect of 
environment preparation for the use of excitonic systems in 
quantum information processing. 
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